SPAGINS: Semiempirical Parameterization for fragments in Gamma Induced Nuclear Spallation* 
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From the empirical phenomena of fragment distributions in nuclear spallation reactions, semiempirical for- 
mulas named SPAGINS were constructed to predict fragment crosssections in high-energy y-induced nuclear 
spallation reactions (PNSR). In constructing the SPAGINS formulas, theoretical models, including the TALYS 
toolkit, SPACS, and Rudstam formulas, were employed to study the general phenomenon of fragment distribu- 
tions in PNSR with incident energies ranging from 100 to 1000 MeV. Considering the primary characteristics of 
PNSR, the SPAGINS formulas modify the EPAX and SPACS formulas and efficiently reproduce the measured 
data. The SPAGINS formulas provide a new and effective tool for predicting fragment production in PNSR. 
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I. INTRODUCTION 


The nuclear spallation reaction is a violent reaction induced 
by light-charged particles (LCPs), such as p,d,t, He, a, 
and uncharged particles, such as neutrons and high-energy 
y. In nuclear spallation reactions, various residual fragments 
are produced. Fragment production cross sections in spal- 
lation reactions are important because they provide consid- 
erable information about the evolving reaction system [1- 
3]. y-induced nuclear reactions, called PhotoNuclear Reac- 
tions (PNRs), have been extensively studied in recent decades 
[4]. Because only electromagnetic interaction occurs between 
photons and nuclei [5—7], PNRs provide a unique tool for 
studying the properties of nuclear forces [8], nuclear struc- 
tural parameters, nuclear astrophysics, and other fields [9- 
11]. The past few decades have witnessed a renaissance 
of experimental PNRs in laboratories, partly because of the 
emergence of new accelerator-based gamma sources that cre- 
ate quasi-monoenergetic photon spectra with high credibil- 
ity based on laser Compton backscattering (LCB) [12]. The 
concept of producing high-energy photons from light photons 
colliding with extremely relativistic electrons was proposed 
by Milburn [13] and Arutyunian et al. in 1963 [14]. Be- 
mporad et al. used a ruby laser scattered with 6 GeV elec- 
trons to produce 425 MeV y photons two years later [15]. In 
2013, the LEPS-II (the laser electron photon experiment at 
SPring-8) produced 1.4 to 2.4 or 1.4 to 2.9 GeV y beamlines. 
With the development of SLEGS (Shanghai Laser Electron 
Gamma Source (SLEGS) at the SSRF (Shanghai Synchrotron 
Radiation Facility (SSRF) [16, 17] and completion of the fu- 
ture SHINE (Shanghai HIgh repetition rate XFEL aNd Ex- 
treme light facility) construction[18, 19], the gamma energy 
of Compton-scattered light sources based on LINAC (LINear 
ACcelator) accelerators will be on the order of GeV. An area 
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of y energy larger than 140 MeV is mainly related to the fields 
of hadron physics, where mesons can be produced [20, 21]. 
GeV photon beamline can be used to investigate subatomic 
and nuclear physics. 

High-energy y rays can induce nuclear spallation reactions, 
which are called PhotoNuclear Spallation reactions (PNSR). 
In recent years, considerable research has been conducted 
on the photonuclear spallation of heavy nuclei [22—24]. In 
1986, the experimental study of PNSR by Shibata et al. [4] 
included the measurement of the cross sections of 24 nu- 
clides in the y+”*‘Cu reaction. The y rays were made from 
bremsstrahlung with maximum end-point energies from 100 
MeV to 1 GeV by counting the irradiated targets, assuming 
a 1/E dependence of the bremsstrahlung spectrum and giv- 
ing a crosssection with a unit of mb [milli-barn/equivalent 
quantum (mb], and the data were used to monitor the flux 
of bremsstrahlung quanta. The cross sections of the frag- 
ments produced in nuclear spallation shared similar character- 
istics, although they were induced by different incident par- 
ticles. Jonsson et al. compared the fragment yields in y and 
LCPs-induced nuclear spallation reactions, and some photon- 
induced spallation cross sections in !?’I were also estimated. 
They conducted experiments on y-induced nuclear reactions 
above | GeV and discussed the principles of high-energy re- 
actions [25-27]. 

The mechanism of spallation reactions was treated as a 
two-step cascade—evaporation process according to Serber et 
al., in which the first step describes the cascade process and 
the second describes the evaporation process [27-29]. In 
the cascade process, the incoming projectile initiates a chain 
cascade by interacting with nucleons inside the target nu- 
cleus, where numerous particles are ejected from the nucleus, 
and residual target nuclei form highly excited hot fragments. 
Then, in the evaporation stage, hot fragments are deexcited 
by evaporating nucleons or nuclear clusters to form the final 
products. However, in the evaporation process, the memory 
of cascade residual formation is lost, which leads to a very 
similar deexcitation process between the photon- and proton- 
induced reactions [30-32]. 

The main improvement in fragment production in PNRs 
was the development of the TALYS toolkit [33]. The TALYS 


toolkit provides a complete and accurate simulation of the 
nuclear reactions of light-incident particles with energies up 
to 200 MeV using an optical reaction model. However, a 
notable difference was found between the experimental and 
TALYS predicted fragment cross sections [34]. The cas- 
cade evaporation model shows that the fragment cross sec- 
tions in PNSR and proton-induced spallation reactions share 
similarities [29, 35, 36], which motivated us to compare 
their fragment distributions. Considering the lack of system- 
atic prediction models for fragment cross sections in PNSR, 
we propose a semiempirical PArameterization for fragments 
in gamma-induced nuclear Spallation (SPAGINS) based on 
available measured data, as well as theoretical guidance from 
the EPAX, SPACS, and TALYS models. 

The remainder of this paper is organized as follows. In 
Sect. II, a brief introduction to the SPAGINS formalism is 
provided. In Sect. III, the construction process of SPAGINS 
formulas and predicted results by SPAGINS formalism are 
compared to Rudstam-fitting data and TALYS predictions, as 
well as the measured data. A summary of the study is pro- 
vided in Sect. IV. 


Il. METHODS 
A. TALYS toolkit 


The latest version, TALYS-1.96 [33, 37], was adopted in 
this study to guide the construction and verification of SPA- 
GINS formulas. The optical reaction model governs the basic 
concepts in the TALYS toolkit[38]. TALYS toolkit includes 
nuclear reactions induced by n,p,d,t, 23He, a, y, and with 
incident energies ranging from | keV to approximately 200 
MeV [39]. The valid range of target nuclei was between 12 
and 339 (12 < A < 339). In TALYS, the nuclear reaction 
process is divided into three stages: (1) the independent par- 
ticle stage, in which the incident particle is partially scattered 
and partially absorbed, similar to light waves passing through 
a translucent glass sphere; (2) the compound system stage, 
which follows the time of the incident particle being absorbed 
by the target nucleus, energy exchange occurs between them, 
and a compound system is formed; and (3) the end stage, 
in which the compound system is decomposed into outgoing 
particles and residual fragments. The calculation of the cross 
section of a specific residual fragment is based on the manual 
of the TALYS toolkit [40]. The default input parameters were 
adopted in TALYS calculations. 


B. Emprical EPAX and SPACS formulas 


The cross sections of the fragment productions with the 
mass number A and charge number Z for spallation reactions 
can be divided into three terms: the total reaction cross sec- 
tion, mass yield distribution, and charge distribution, which 
were first proposed by Rudstam [41, 42], the five-parameter 
fitting formula, and the development of empirical EPAX and 
SPACS formulas. The five parameters in the five-parameter 


fitting formula (named the Rudstam formula in this arti- 
cle) were obtained by fitting with the nonlinear least-squares 
method; thus, it cannot be used to predict nuclides. 


The main characteristics of fragment production in spalla- 
tion reactions were included in the EPAX parameterization (a 
universal empirical parameterization of fragmentation cross- 
sections), which was proposed by Siimmerer et al. in 1990. 
The EPAX formulas inherit the ideas of Rudstam and Sil- 
berberg [43], who aimed to describe the fragmentation of 
medium-to-heavy mass projectiles. The fragment cross sec- 
tion was independent of the incident energy of the reaction 
system above 140 MeV/u. The updated versions of EPAX2 
and EPAX3 successfully describe the production of frag- 
ments in projectile fragmentation reactions above 100 MeV/u 
[44, 45]. The improved formulas for fragments produced in 
nuclear spallation reactions based on EPAX have been pro- 
posed by Schmitt et al. in 2014, which is named SPACS (a 
semiempirical parameterization for isotopic spallation cross 
sections) [46, 47]. In the SPACS formulas, more than 40 
parameters were adopted to reestablish the mass yield and 
explicitly state the dependence of the yields of the residual 
fragments on the bombardment energy [48]. The charge dis- 
tribution of the fragments was given by EPAX [44, 45, 49]. 
The dependence of the fragment cross section on the collision 
energy, shell structure, and even-odd effect was further con- 
sidered in the SPACS formulas. The readers can refer to Refs. 
[44, 45] for detailed descriptions of fragments cross sections 
in EPAX and to Ref. [47] for SPACS parameterizations. The 
main formulas in both EPAX and SPACS were adopted to 
construct the SPAGINS formulas in this study, which are in- 
troduced in Sect. II C. 


C. Phenomenological Isotopic Distributions in PNSR 


In this subsection, the fragment isotopic cross sectional dis- 
tributions are compared to determine the basic ideas for devel- 
oping the SPAGINS formulas. In the first step, we illustrated 
the similarity in fragment production between ~y and proton- 
induced nuclear spallation reactions. The TALYS-1.96 was 
adopted to predict the isotopic distributions for the y+°*Cu 
and p+°?Ni reactions at 100 and 200 MeV, where the mass 
and charge numbers of the reaction systems were the same. 
In Fig. 1(a) and (b), the isotopic cross section distributions 
in the y+°%Cu reaction (open symbols) share the same pat- 
tern as those in the p+°?Ni reaction except they have smaller 
magnitudes. The quantity is defined as the ratio of the frag- 
ment cross section in proton-induced spallation (a(f)») to y 
induces spallation (o(f)1), Rp» = o(f)p/a(f)y- The Rp + 
values for the isotopic crosssections in p+®Ni and y+°Cu 
reactions are shown in Fig. 1(a) and (b), respectively. The 
Ry /-, values for isotopes in the reactions at an incident energy 
of 100 MeV were approximately 150 and 400 at an incident 
energy of 200 MeV. This indicates that the formula for frag- 
ment cross sections in the proton-induced nuclear spallation 
reaction can be adopted for y-induced ones by further consid- 
ering the incident energy dependence. 
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Fig. 1. (Color online) The calculated isotopic cross section for fragments by TALYS-1.96 in the y+°?Cu and p+°Ni reactions at E„ = 100 
MeV [in (a)] and 200 MeV [in (b)], respectively. The ratio of the isotopic cross sections in the pt? Ni and yt Cu reactions are plotted in 


the inserted figure. 


D. Rudstam Formula 


The Rudstam formula has two types of distributions for 
determining the fragment crosssection in PNSR. One uses 
charge distribution and mass yield distribution (CDMD), and 
the other uses isotopic distribution and elemental distribution 
(IDED)[41]. In this study, the CDMD in Refs. [4] is adopted, 
and the Rudstam formula reads, 


6PW?2/3 


0(4, 4) = T79(ePA T) 


exp(PA—W|Z—SA+T A?|9/2), 
(1) 


where P, W, S,T and ô are free parameters. P defines the 
slope of the mass yield curve, and W defines the charge dis- 
tribution width. S and T describe the most probable charge 
and define the peak locations of the charge distribution (or 
isotopic distribution), respectively. ¢ denotes the total inelas- 
tic yield of the reaction. The parameters were determined by 
performing a nonlinear least-squares fit to the measured frag- 
ments in y+""'Cu PNSR within the range of Æ, from 100 to 
1000 MeV [4] and adopted in this study to investigate the in- 
cident energy dependence of the fragment cross section (see 
Sect. III A). 


In the second step, with the help of the Rudstam formula, 
the fragment cross sections of chromium isotopes were com- 
pared for the y+°Cu reaction at incident energies ranging 
from Æ, = 100 to 1000 MeV. 


HI. SPAGINS FORMULAS 


In this section, the procedure for developing the SPAGINS 
formulas is first described, and then the SPAGINS predictions 
are verified to fragment the cross sections in PNSR through a 
comparison with the measured results. 


A. Developing SPAGINS Formulas 


For a residual fragment (Z, A) produced in fragmentation 
or spallation reactions, the cross section can be described as 


o(A, Z) = orY(A)Y (Zprob — Z) a, (2) 


where Op is the normalized total reaction cross section [50, 
51], and Y(A) is the mass yield. The third term represents 
the isobaric yield for a given mass A. The construction of 
SPAGINS formulas started with this formula to describe frag- 
ment production. In this study, modifications were made to 
the mass and charge distribution terms. 

The terms that influence op are the masses of the projec- 
tile Aproj and target Atar. The energy-dependent function ôg 
considers the effects of transparency, Pauli blocking, and the 
Coulomb barrier B. The energy of the collision system is in 
the center-of-mass framework, and the quantity Xm corrects 
the intensity of the optical model interaction at low energies 
[46, 47]. Because photonuclear reactions involve only elec- 
tromagnetic interactions, the Coulomb barrier is zero. The 
mathematical expressions for these quantities are given in the 
SPACS formalism [5, 6]. 
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Fig. 2. (Color online) The Rudstam fitting formula according to Eq. 
(1) for production cross sections of Cr isotopes in y+°Cu reactions 
from 100 to 1000 MeV. 


The measured fragment cross sections in the y-++Cu spalla- 
tion reaction [4] were used as an example to perform the anal- 
ysis. Natural copper has two stable isotopes: 69.17% °3Cu 
and 30.83% °°Cu. The findings of this study revealed that the 
isotopic effects in the fragment cross sections were minimal. 
Thus, °*Cu was adopted as the spallation target instead of nat- 
ural copper, which is mainly composed of ĉ3Cu and ®°Cu. To 
maintain the same masses and charge numbers in spallation 
systems, a proton-induced system of p+°Ni was selected for 
comparison. 


The calculated fragment cross sections were compared to 
the measured cross sections according to two main aspects. 
The first method considers the incident energy dependence of 
the fragments. With an increase in energy, the cross sections 
of the residual fragments produced by PNSR increase expo- 
nentially. The second factor is the charge number Z. For the 
generation of different types of residual fragments, the differ- 
ence between the calculated residual fragment cross section 
and experimental value constantly changes with the change 
in charge number. To reduce the difference between the sec- 
tion value calculated using the empirical formula and the ex- 
perimental value, the corrections for energy and charge are 
necessary. 

Schmitt et al. found that the mass distributions of the spal- 
lation and fragmentation reactions cannot be described by 
the same mathematical expression [46]. Importantly, in spal- 
lation reactions, the energy dependence of Y(A) cannot be 
ignored[44, 45, 48]. Schmitt et al. added an energy depen- 
dence related to Y (A) to the SPACS formalism [43, 45, 49]. 
The energy dependence of Y (A) is discussed in Sect. IIIB. 

The third term in Eq. (2), Y(Zprop — Z)|4 describes 
the isotopic distribution, which is insensitive to the interac- 
tion mechanism but is mainly controlled by the level den- 
sity. In this study, the application ranges for certain quantities 
were modified. The magnitude of the brute-force factor was 
changed, and the charge dependence was modified, as shown 


in Sect. III B. 


B. Main Formulas in SPAGINS 


Mass yield Y(A) in Eq. (2) is divided into two parts con- 
sidering the fragment contributions of the central collisions 
Y(A)cent and peripheral collisions Y(A)prph, which are ex- 
pressed as 


Y(A) = Y (A) cent ale Y (A)prph- (3) 


Y (A) depends exponentially on the incident energy. Particu- 
larly, the production of Cr isotopes was studied to determine 
their dependence on the incident energy. To observe how the 
fragment cross section depends on E4, the Rudstam formula 
in Eq. (1) with the parameters according to Ref. [4] is adopted 
to estimate the Cr isotopic productions in y+°*Cu reactions 
with Æ, from 100 MeV to 1 GeV (see Fig. 2). Within the 
range of 100 MeV< E, <500 MeV, the isotopic distribution 
increases with the incident energy, whereas it becomes very 
similar when Æ, > 500 MeV. Thus, the trend of Y (A) is pa- 
rameterized as the following equations for different incident 
energy ranges: 


Y'(A)=G-Y(A). z 
G depends on the incident energy in the form of 
10(27—100)/100 100 < E, < 220; 
220 < EY < 500; (5) 
10(#7—500)/2000 500 < E, < 1000. 


G = } 19(27~-220)/300 


E, is sorted into three different ranges: (1) from 100 to 
220 MeV, (2) from 220 to 500 MeV, and (3) from 500 to 1000 
MeV, which empirically reflect the energy dependence of the 
isotopic distributions shown in Fig. 2. 

In the EPAX formulas, Zprob describes the deviation of the 
most probable charge from the position of G—stability valley 
(Zg) by the quantity A. For photonuclear spallation reaction, 
the parameters A used in EPAX are as follows: 


A= 454°, (6) 
A’ = All + dı (A/Aproj — d2)°]. (7) 
The SPACS formula modifies the parameter R in EPAX, 
which is also the formula used: 
R= Rj? RP” 
Ro? = To exp[r4(Zproj = ZBp)l, (8) 
RY = exp{— In| R1 (A/t2)]}. 


The residual fragments in the spallation reaction can be di- 
vided into proton- and neutron-rich fragments, which are pa- 
rameterized in different forms. The SPACS refined the nor- 
malization of the charge-dispersion curve based on EPAX3, 
and the so-called “brute-force factor” fa (fp) is used to adjust 
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Fig. 3. (Color online) The isotopic cross section in the 100 MeV (a), 130 MeV (b), and 160 MeV (c) y+%*Cu reactions. The SPAGINS 
predictions are plotted as open symbols, and full lines represent the fitting lines by the Rudstam formula, which are used to guide the eye. 
The measured fragments in the y+"“*Cu reactions (taken from Ref. [4]) are plotted in full symbols. 


the distribution of neutron-rich (proton-rich) fragments. For 
(Zs — Z) > (Zproj — Zap + b2), 


fa = 10091] Zw Zap (Zp—Z+ Zone Zap tb2)°) (9) 
Otherwise, fa = 1. For Z > Zexp 
h= 1/[109£/42](2- Ze»), (10) 


When the charge of the residual fragment was less than that 
of the target nucleus, the SPACS predictions agreed with the 
measured data, whereas when the charge of the residual frag- 
ment was close to that of the target nucleus, the SPACS pre- 
dictions differed significantly from the measured values. In 
addition, fn and fp determine the peak position of the Gaus- 
sian distribution of the isotopes. The following conditions of 
“brute force factors” are formulated for those residual frag- 
ments on the valley of @ stability. When Zg > Z — 0.75, fn 
is used, and fp is used in all other cases. For these reasons, 
in y+°4Cu spallation reactions, the residual fragments with 
charge numbers ranging from 21 to 29 are divided into three 
groups. The first group consisted of residual fragments with 
Z < 21, second group consisted of those with 21 < Z < 26, 
and last group consisted of those with Z > 26. The parame- 
terized brute force factor is described in Sect. VI. 


C. Validation for SPAGINS formulas 


In Fig. 3, the predicted isotopic cross section distributions 
by the SPAGINS formulas in y+°*Cu reactions and solid line 
used to guide the eye is the five-parameter fitting formula 
from Rudstam. The cross sections for fragments in the corre- 
sponding y+""Cu results are plotted for comparison in solid 
symbols, for which the incident energies Æ, range from 100 
MeV to | GeV at 100, 130, 160, 220, 310, 400, 500, 800, 
and 1000 MeV. The mass (charge) number of the fragments 
ranged from 38 (19) to 64 (29). Owing to the similarities in 
the types of graphs, only Æ, values at 100, 130, and 160 are 


displayed in this study. The SPAGINS formula reproduced 
the measured data well, except for those near the target iso- 
topes. 


The excitation function of the residual fragment, that is, the 
dependence of the fragment crosssection on the incident en- 
ergy of the reaction, reflects how the probability changes with 
E. Figure 4 shows the cross sections of fragments from 2K 
to ĉtCu, which are produced in the y+°*Cu reaction within 
E. from 100 MeV to 1 GeV, in which both the SPAGINS 
predictions and Rudstam fitting data are compared to the mea- 
sured results of the y+"Cu reaction. Both the SPAGINS and 
Rudstam formulas can reproduce the experimental excitation 
functions in y+""'Cu reactions, except that the SPARGINS 
prediction underestimates the measured function for °’Ni. In 
the SPAGINS formula, Eq. (26) in the APPENDIX, and the 
brute force factors fn and fp the fitting data of the Rudstam 
formula are plotted in depend on E, and influence the frag- 
ment cross section [46] when E, is changed. In the Rudstam 
formula, all five parameters depend on incident energy [41]. 
The good agreement between the SPAGINS predictions and 
measured results for the excitation curves of fragments sug- 
gests that the parameterization of the incident energy depen- 
dence of fragment crosssections reflects the inner mechanism 
of fragment production in PNSR. The mean relative errors are 
shown in Figs. 5. The dashed line indicates that Gexp/o¢al 18 
equal to 1. 


To further verify the SPAGINS formulas, the predicted 
cross sections for fragments using the TALYS toolkit and 
SPAGINS formulas in the 100 MeV y+°*Cu reaction are 
compared in Fig. 6. The mean relative error between the 
predictions of the SPAGINS formula and experimental value 
in the 200 and 900 MeV y+°°Co reactions [52] is shown in 
Figure 7. In general, the data predicted by both TALYS and 
SPAGINS are consistent with the measured data, except for 
the neutron-rich fragments, in which TALYS predicts rela- 
tively lower data. If this phenomenon is valid for neutron-rich 
fragments, then the SPAGINS predictions agree with the mea- 
sured data in PSNR. Based on comparisons between the SPA- 
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Fig. 4. (Color online) Excitation functions for fragments 22K, Hsc, y, 51Cr, 52Mn, 59Fe, 56C0, 57Ni, and Cu produced in reactions of 
y+*Cu with E. from 100 MeV to 1 GeV. The measured data in y+""Cu reactions (solid circles) are taken from [4]. The predictions by 
SPAGINS and the fitting data of the Rudstam formula are plotted in open squares and open triangles, respectively. 
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Fig. 5. (Color online) Mean relative error of 42,43 K 44,46~48 Sc, 484, 49,51 Cr, 52:54:56Mn, 59Fe, 56~58:60C0, STNG and 8%61:64Cu produced 
in the reactions of y+6*Cu for which the incident energies E, ranged from 100 MeV to 1 GeV. The dashed lines denote Gexp/Gcai is 1. 


GINS, Rudstam formula, and TALYS toolkits, the SPAGINS IV. SUMMARY 
formulas provided reasonable predictions for fragments pro- 


Cee an PANK, Considering the lack of an effective model to predict 


the fragments produced in PNSR for high-energy y rays, 
semiempirical formulas named SPAGINS have been pro- 
posed, which are suitable for PNSR within the range of 100 
MeV < E, < 1 GeV. The following procedure was fol- 
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Fig. 6. (Color online) Production cross sections for Cr, Mn, Fe, Co, Ni, and Cu isotopes in the 100 MeV yt+%Cu reaction. The measured 
data for y+"'Cu reactions (solid circles) are taken from Ref. [4]. The predicted results by the TALYS model and SPAGINS formulas are in 


open triangles and open squares, respectively. 


lowed to construct the SPAGINS formulas: (1) TALYS-1.96 
was adopted to find the similarity of fragment production 
in photon- and proton-induced spallation reactions with the 
same mass and charge numbers, that is, y+°%Cu and p+°?Ni 
at 100 and 200 MeV. In this step, the isotopic cross sectional 
distributions of the different elements had similar shapes but 
different magnitudes. This enabled the borrowing of the 
main concepts of the models for proton-induced spallation 
reactions. (2) The isotopic cross sectional distribution of 
chromium in y+°4Cu at E. from 100 MeV to 1 GeV fit- 
ted from the Rudstam formula was used to parameterize the 
incident energy dependence of mass yields in PNSR reac- 
tions. (3) Based on steps (1) and (2), the SPAGINS formu- 
las were constructed based on the SPACS formulas for light- 
charged-particle-induced nuclear spallation reactions, as well 
as EPAX formulas for projectile fragmentation reactions, by 
implanting proper modifications to describe fragment distri- 
butions. The main characteristics of the excitation function 
of the mass yield were also established in this step. (4) The 
predictions of isotopic cross sections by SPAGINS formulas 
were compared to the measured data and Rudstam formula, 
which were in good agreement. The excitation curves of the 


fragments also support the conclusion that the SPAGINS for- 
mulas can predict the fragment cross section at different E, 
values from 100 MeV to 1 GeV. 


Compared to the Rudstam formula, the SPAGINS formu- 
las overcome this limitation because the parameters are de- 
termined from a series of measured data and restricted to lim- 
ited reaction systems. Meanwhile, the SPAGINS formulas 
also expand the applicable range of the incident energy com- 
pared with the TALYS model (lower than 200 MeV) to 100- 
1000 MeV. The success of SPACS and EPAX formulas in de- 
scribing fragment distributions makes the SPAGINS formulas 
yield reasonable fragment cross sections, and the fragment 
excitation curves are suitable for PNSR from E, from 100 
MeV to | GeV. Currently, SPAGINS calculations range from 
Fe to Zn targets, with charge numbers between 26 and 30 
and mass numbers between 58 and 68, which cover common 
metallic materials in nuclear industrial applications. Consid- 
ering the rapid development of high-energy y-ray facilities, 
SPAGINS formulas provide an effective method for estimat- 
ing fragment production, y-nuclear activation, and radiation 
protection. 


Fig. 7. (Color online) Mean relative error of 98598599 Cy. REK, 
44,46~48 Sc, 48~51. Cr, 52:54:56Mn and 52:53 Re produced in the 200 
MeV and 900 MeV 7+°°Co reactions. The measured data are taken 
from [52]. The dashed line means that Gexp/Ocat is 1. 


V. APPENDIX I: SPAGINS FORMULAS 


Many SPACS and EPAX formulas are incorporated into the 
SPAGINS formulas. Owing to the definition of the individual 
parameters in the formula, SPAGINS adopts gamma quanta 
as the target. To present the SPAGINS formula more clearly, 
detailed formulas are provided in this section. 

For a fragment with mass and charge numbers (A, Z), the 
production cross section in PNSR is 


o(A, Z) = oRY (A)Y (Zprov — Z) A. (11) 
For the first term, the reaction cross section is 


or = 10mr2(A2,. + AŠ, + bp)2(1— B/Eom)Xm, (12) 


proj 
where rp = 1.1 fm, and 

B = 1442 p10j Zar H, (13) 

Eom = E proj Aproj / (Aproj + Atar), (14) 


1 1 A 
H = fT proj + Ttar + 1.2(AS, + Ags) / Ean. (15) 


proj 


One can determine fproj = 1.29r proj» and Ttor = 
3 -5 -5 rms __ 
0.891A5;(1 + 5.565 A oj — 1.044 gi) Tar = 0.85. 
i 
dg = 1.85S + (0.16S9/Ečn) — Cp, (16) 


where 


S = A3,A3, /(A3, + A3 (17) 


proj proj) 
and 
Cp = 2.05[1 — exp(—Eproj/40) — 0.292 exp(— Eproj/792) 
cos(0.229E2.453)] (18) 
For Aproj < 200, 
Xm = 1 — xı exp|— Eproj/ (Xıs1)], (19) 


in which, 


xı = 2.83 — 3.1 x 107? Aproj + 1.7 x 1074A? (20) 


proj? 


and 
sı = 0.6, for Aproj <12; (21) 
s1=1.6, for Apo = 12, (22) 
sı = 1.0, for Aproj > 12. (23) 


For the second term, the mass distribution considers the con- 
tributions from the central (Y(A)cent) and peripheral colli- 
sions (Y (A)prph), 


Y(A) = Y(A)cent F Y (A)prph- (24) 


When considering the incident energy dependence of the 
mass yield, 


Y'(A)=G-Y(A) (25) 
where G is sorted into the following three regions: 


10(#7=100)/100. 100 < FE, < 220; 


220 < Ey < 500; (26) 
10(21—500)/2000 599 < E, < 1000. 


G = 2 19(£1—220)/300 


The contribution of central collisions to a specific fragment 
cross section is 


Y (A) cent = (Atength/{1 le exp|( Aproj = Aza 
—A) / Acent-fluct] })/ Avent: (27) 


where 


Acent = cent/ €XP(Aproj/Feent); 
Avent = Aproj {1 — exp[— In(A cent) (Eproj /1000)*"]}, 
Atength = 1, 
Acent-fuct = atut AGE ( Ero / 1000), 
Afluct = Qfluct + Bact Aproj- 


The contribution of peripheral collisions to a specific frag- 
ment crosssection is 


Y(A)periph = Aperipn exp{[A — (Aproj—Aair)]/Aaitt} 


x (Eproj / 1000) Pret Aii (28) 


where 
Aperiph = (periph F Bperiph Aproj 
Aagitt = Oise + Baitt Aproj 
proj 


The third term o (A, Z) is parameterized by the number of 
neutron- and proton-rich fragments. For neutron-rich ones, 


Y (Zprob z= Zla = nfn exp(— R| Zprob Se Z|""), (29) 
and proton-rich 
Y (Zprob — Z)|a = Nfp exp(—R|Zprop — Z|“”). (30) 
Zprob and Zg take the form of 
Zprob = Zg + AGP +0.002A+A’, (31) 
_ A 
ZB = T98F0.0155A773" (32) 
For a neutron-rich projectile, 
Am = {n1 (A/Aproj) + 2[(Aproj — A) /Aproj]}” 
x (Zproj = Zep); (33) 
For proton-rich projectiles, 
Ah = [exp(p1 + p2A4/Aproj)](Zproj — Z6p) (3.4) 


in which A = AsA?, A’ = A[1 + di(A/Aproj — d2)?], n = 


y R/T, and R = Ry? RY, 
For neutron-rich spallation targets, 


Ro = ro exp[r3(Zproj — Zap); (35) 
For a proton-rich spallation target, 

RE = ro exp|ra(Zproj — Zap)]- (36) 
One also has, 

RP = exp{—In[Ri(A/tz)]}. (37) 


For neutron-rich fragments, 
U” = Un + Un A/Aproj; (38) 
and proton-rich fragments. 
UP = Unt + Up2 Aproj- (39) 


The brute-force factors f,, presented in Sect. VI. 


VI. APPENDIX II: BRUTE FORCE FACTOR 


The “brute force factor” in Eqs. (9) and (10) were param- 
eterized according to the incident energy of the reaction and 
different ranges of fragment charge numbers. For the energy 
range of 100 MeV < E, < 220 MeV, when Z < 26, 


_ 1 
h= 10(9F/42)(Z— Zexp) 
x 1 910-5(29- Zpo)][5.575—0.175Z] 


x 1gl-12+0.5(4-1)] 


(40) 
fa = 104 x jot 23+ (4-0) 
x 100-529- Zpro;)][4.7—0.12] 


where 
dF/dZ = 1.2 + 0.647(4/2)®™?, (41) 
Zexp = Zprob + AF'/dZ In(10)/(2R); 
when Z > 26, 
f= 1 x 10 24+(Z-1)] 
P 1 Q(4F/dZ)(Z— Ze) , (42) 


fa = 104 x 100-425+(2-) 


For the incident energy range of 220 MeV < E, < 500 MeV, 
when the charge of the fragments Z < 21, 


1 


[5.5—0.3(Z—1)] 
fo = T00F/42Z Za) * 10 


" (43) 
A = 104 x 1gl-15.5+0.75(Z—1)] 


i 


when 21 < Z < 26, 


1 
= x 10178:2+0-4(7-1)]. 
fo 1Q(dF/dZ)(Z—Zexp) (44) 


fa = 104 x 101- +2-6+0.6(7-1)]. 


when Z > 26, 


1 


_ [-10+0.5(Z—1)] 
f= 10F/AZ Z-Z) * 10 


a 


(45) 
fa = 104 x 1027-7+1.2(2-1)], 

For the incident energy range of 500 MeV < E, < 1000 
MeV, when Z < 21, 


1 
= [5.7—0.3(Z—1)] 
fo = OFAI ZZ) * 10 


x 1 10-5(29— Zprej)][2-2-0.12] 


(46) 
fa p 104 x 10[719-6+0.6(7—1)] 
x 1910-5(29—Zproj)][1-15—0.05Z]. 
when 21 < Z < 26, 
1 
= [—7.4+0.4(Z—1)] 
h= 10F/aZ Z-Z) ” 10 
0.5(29—Zproj)]|[-2+0.1Z 
x 1910-56 MI l (47) 
fa = 104 x 10 9-5+0-5(2-1)] 
x 1910-5(29—Zorej)][-1-9+0.1Z] 
and when Z > 26, 
f = 1 x 10[79-8+0-5(2-1)] 
P 10(4F/d42Z)(Z— Zexp) : (48) 


fa = 104 x 10[714-2+0.7(7-1)], 


where parameters from Eq. (3) to Eq. (48) can be found in 
Refs. [46, 47]. 
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